On the Stability of the Compacton Solutions 
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Abstract 

The stability of the recently discovered compacton solutions is studied by 
means of both linear stability analysis as well as Lyapunov stability criteria. 



From the results obtained it follows that, unlike solitons, all the allowed com- 
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' pacton solutions are stable, as the stability condition is satisfied for arbitrary 

values of the nonlinearity parameter. The results are shown to be true even 
for the higher order nonlinear dispersion equations for compactons. Some new 
conservation laws for the higher order nonlinear dispersion equations are also 
presented. 
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The observed stationary and dynamical patterns in nature are usually finite in extent. 
However, most of the weakly nonlinear and linear dispersion equations so far studied admit 
solitary waves that are infinite in extent, although of a localised in nature. Therefore, 
the recently discovered compacton solutions (i.e. solitary waves with compact support), of 
the nonlinear dispersive K(m, n) equations have become very important from the point of 
study of the effect of nonlinear dispersion on pattren formation as well as the formation 
of nonlinear structures like liquid drops etc. The compacton speed depends on its height, 
but unlike the solitons, its width is independent of the speed, a fact, which seems to play a 
very crucial role in its stability property. Compactons have remarkable soliton like property 
that they collide elastically. However, unlike soliton collisions in an integrable systems, the 
point where two compactons collide is marked by the creation of low amplitude compacton- 
anticompacton pairs [1,2]. In fact, it is now known that the K(m,n) system of equations 
are not integrable [1,2]. This suggest that the observed almost elastic collisions of the 
compactons is probably not due to the integrability and thus the mechanism responsible for 
the coherence and robustness of the compactons remains a mystery. Stability analysis of the 
compacton solutions may provide some clues in this direction. As has been said above, the 
stability of the compactons is crucial in the context of its applications in the study of pattern 
formation. Beside, the stability problem of the K(m, n) equations is interesting because, for 
such equations with higher power of nonlinearity and nonlinear dispersion, the phenomena 
of collapse is possible. Also, in the context of soliton equations, from the stability analysis it 
has been shown that the higher order linear dispersion term stabilises the solitons [3,4]. In 
this regard it will be of interest to see what role does the higher order nonlinear dispersion 
term plays with respect to the stability of the compacton solutions of theK(m, n) type 
equations. 

In this communication we report on the stability analysis of the compacton solutions of 
the nonlinear dispersion K(m,n) type equations as considered by Cooper et al [2]. We use 
both the linear stability analysis and the Lyapunov stability criteria to analyse the problem. 
We start with the K(l,p) equations 
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u t + u x u l 2 + a[2u 3x u p + Apu p l u x u 2x + p(p - l)u p 2 u x ] = (1) 

These equations have the same terms as the K(m, n) equations considered by Rosenau et 
al [1] but the relative weights of the terms are different, leading to the fact that, whereas 
the K(l,p) equation (Eqn.(f)) can be derived from a Lagrangian, the K(m,n) equation 
considered in [1] donot have a Lagrangian. For the sake of comparison, it may be noted that 
the set of parameters (m,n) in [1] corresponds to the set (I — l,p+l) in Eqn.(l) [2]. Assuming 
a solution to Eqn.(l) in the form of a travelling wave u(x,t) = m(£), where £ = x — Dt, 
Eqn.(l) reduces to the same K(l,p) equations considered by Cooper et al (Eqn.(7) in [2]) 
for the compactons solutions. The compacton solutions to Eqn.(l) are given by [2] 

«(0 = [g(p+l)(p + 2)]W[ } (2) 

2 2Ja(p + l)(p + 2) 



in case I = p + 2, < p < 2 and for | £ |< ^ya(p + l)(p + 2) and tt(£) is zero otherwise. 
Note the width of the compacton is independent of its speed (amplitude). The fact that the 
width of the compacton will always be independent of its amplitude in case I = p + 2 follows 
from the invariance of Eqn.(l) under the scaling transformation x — > ax,t — > [3x,u — > jx. 
We now consider the stability of these compacton solutions. 

Linear stability- Eqn.(l) can be obtained from the variational principle 5(H + DP) = 0. 
Using the relations 

/+oo f+OO 
u n dx and J 2 — u p u 2 x dx (3) 
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we can write the corresponding conserved Hamiltonian and momentum as 

(r+tto+ty Pc=l 2 h (4) 

Using Eqn.(l) and the equation obtained from the scaling transformation x — > /3x, we get 

/ p+2 _( i + P )DP c md aj2 = J^ (5) 



(p+l)(p + 2) 2(p + 2) 2(p+2) 
so that H c = —2DP c /(p + 2). Now considering the general scaling transformation u 
/i 1 / 2 -u(Aa;), H c and P c are transformed to H(X,fj) and P(X,fj) and we get 



P+2 

H x,, ) = a x^j 2 - 1 j-tl— Ip+2 + E D p c (6) 

where $(A, //) = H(X,/jl) + DP(X, /i). The equations|| = |^ = gives the stationary point 
at A = = 1 (compacton equation) and near this point, using the Taylor series we get for 
fj, — A ( the transformation in that case does not change the momentum P) 

»(A) = SH W = (A - m a(P + 2)i i + 4> ^ - P<P P ' l (7) 

which has a definite sign. If it is positive (negative) the expression 

has a minimum (maximum) at A = 1. 

Now, let us assume that u — u c + v, where | v \« 1 and the scalar product (u c , v) = 0. 
Substituting this in Eqn.(l) we get after linearization 

8tv = d^Lv (9) 

where£ = x — Dt and T=t and the operator L is given by 

L—[D — u 1 - 2 - 2apu p ~ 1 u 2 t - 2avPd\ - ap(p - l)u {p - 2) u\ - 2apu {p ~ 1) u i .d i .} (10) 

One can now run through the steps as given by Karpman [3] (see his Eqns. (29) to (32)) 
and show that as in his case, in our case too the sufficient condition for stability is that the 
scalar product (tp, Lip) > where the operator L is given by Eqn.(10) while ip is a function 
in the subspace orthogonal to u c . However, condition (ip,Lip) > is also associated with 
the extremum of H + DP, since, using the relation 5(H + DP) = one can show that the 
second variation of H{u) and P{u) at u = u c is given by 



5 2 (H + DP) Uc = i / (v,Lv)dC>0 (11) 
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where the operator L is given by Eqn.(10). That is, if the condition (ip,Lip) > is fulfilled, 
then H{u) + DP{u) has a minimum at u = u c . Inversly, the minimum of H{u) + DP{u) at 



u = u c is a sufficient condition of compacton stability with respect to small perturbation. 
Using Eqn.(8) we obtain the condition for the minimum of the perturbed Hamiltonian H(X) 
at A = 1 as p+2 > p — 2 which is obviously true for any p. Thus we see that the condition for 
the Hamiltonian minimum (and hence the sufficient condition for the compacton stability) is 
satisfied for arbitrary values of the nonlinear parameter^. This is unlike the soliton stability 
results, where it has been shown that the stability condition puts a restriction on the allowed 
values of the nonlinear parameter [3,4,6]. Note however that compacton solutions exist only 
for p < 2. 

A la Karpman [3], one can show that the sufficient condition (ip, Lip) > is also equivalent 
to the condition 

(§)>0 (12) 

From Eqns.(2) and (4) one can easily show that the sufficient condition for compacton 
stability (Eqn.(12)) is satisfied for arbitrary values of the nonlinearity parameter p. It should 
be noted that this result is completely in contrast to the usual soliton stability results. It 
is not difficult to see why this is so for the compactons. From Eqn.(2) we see that the 
width of the compacton solutions is independent of its speed (amplitude) D, and the generic 
form of such compactons is u c (£) = AD b cos(c^), where the constants A,b and c depends 
on the nonlinearity parameter. Hence P c = D 2b K, where K is D-independent. Therefore, 
dP/dD > trivially since b > 0. On the other hand, if the width depends on speed (as in 
the case of solitons [3,4]) with the generic form of the solution as w(£) = AD b cos(cD a l;), then 
dPj dD > only if 2b > a which will depend on the particular theory (soliton equations) . It 
should be noted that the above stability condition (Eqn.(12)) is obtained by assuming that 
there is only one -ve eigenvalue for the operator L (Eqn.(10)). The validity of this conjecture 
has been proved from numerical experiments for many other systems, like the third and fifth 
order KdV equations as well as the nonlinear Schrodinger equations [3] . At present we donot 
have any evidence to show that this conjecture is also valid for our operator L (Eqn.(10)), 
except for the fact that the result that follows from using this conjecture also agrees with 



the result as obtained from the Hamiltonian minimum condition (Eqn.(ll)) as well as an 
independent analysis of the stability by Lyapunov method as shown below. 

These results are also true for the higher order dispersion equations. For example, 
consider the fifth order nonlinear dispersion equations, the K(m,n,p) equation [5,7], for the 
compactons. The Hamiltonian for such system is respectively given by [5] 



As has been shown in [5], the compacton solutions corresponding to the Hamiltonian in 
Eqn.(13) are allowed for the nonlinearity parameter k = m = n = p in the range 2 < k < 5 
thereby meaning that the effect of the higher order nonlinear dispersion term is to increase 
the range of the nonlinearity parameter for which the compactons solutions are allowed. 
Considering the small perturbation u = u c + v as before, we can show that [8] even for these 
higher order nonlinear dispersion equations , the sufficient conditions for the compacton 
stability as given by Eqns.(ll) and (12) is satisfied for arbitrary values of the nonlinear 
parameter k. 

Lyapunov stability- The above theory of linear stability analysis for the compactons is based 
on the linearization of equations for compacton perturbations. This method has some in- 
herent limitations connected with the linearization. So we present another approach to 
the stability problem based on the Lyapunov method, which, instead of linearization, uses 
sharp estimates. The effectiveness of this method has been demonostrated by Weinstein [9] 
and others [10]. In this method of analysis, it is sufficient to prove that the Hamiltonian 
is bounded from below for fixed momentum P and the compacton realizes the Hamilto- 
nian minimum. We consider here the stability of the compacton solution of Eqn.(l). From 
Eqn.(3), we have 




(13) 



and [7] 




(14) 




(15) 
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Also, 

max {u^~) < - J | u p / 2 u x || u \ dx (16) 
Using Holder's inequality, we get 

max (u^) < ?±±[J ' vPu 2 x dx] l / 2 [j ' u 2 dx} 1 ' 2 < P -±^ J^ 2 (2P) 1 ' 2 (17) 

From Eqns. (4) we then have 

min l v a. 4 2p _e_ 2 P +4 

* ^ * [<•* " (p+1)(p + 2) (^)^^(2P)^] (18) 

Thus is bounded from below. On calculating the minimum of the right hand side we 
find H min = —^aJ 2 Thus we see that H is bounded from below for arbitrary values of the 
nonlinearity parameter p. Now , from Eqns. (4) in Eqn.(5) we can immediately see that 
H c = H min = — 2 p+2) ' ^- e -' ^ ne com P ac ton realizes the Hamiltonian minimum and hence this 
proves the stability of the compacton solutions in the Lyapunov sense. Thus, the Lyapunov 
stability analysis show that all the allowed comacton solutions (i.e. p < 2) are stable, since 
the condition for boundedness of the Hamiltonian and H c = H min are valid for arbitrary 
values of the nonlinearity parameter p. It can be shown that [8] we also get similar result 
from the Luapunov stability analysis of the higher order nonlinear dispersion equation as 
given by the Hamiltonians in Eqns. (13) and (14). 

Before concluding, we would like to mention that in our earlier paper [5] we had reported 
only one conservation law for the higher order K(m, n,p) equations given by (Eqn.(2) in 



ut + (3\{u m ) x + P 2 (u n ) 3x + P 3 (u p ) 5x = 0, m, n, p > 1 (19) 

We now find that, like the K(m,n) equations as considered by Rosenau et al [1], the higher 
order K(m,n,p) equations also have four conservation laws for m = n = p, with the same 
conserved quantities as for the K(m, n) equations, i.e. 

Qi = u, Q2 = u m+1 , Q3 = ucosx and Q± = usinx (20) 
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We have checked that even for the seventh order nonlinear disperison 

K(m,m,m,m) equation (see Eqn.(53) in [5]) there are four conservation laws as 
above. We suspect that even the generalised arbitrary odd-order nonlinear dispersion 
K(m,m,m,m, ....) equations (see Eqn. (36) in [5]) may also support similar four con- 
servation laws. It should however be noted that for the K(l,p) equations (Eqn.(l) above) 
and its corresponding higher order K(m,n,p) equation (Eqn. (13), (14) above), which are 
derivable from a Lagrangian and whose stability problem is considered here, there are only 
three conservation laws [2,5]. 

To conclude,we would like to point out the important difference between the soliton and 
compacton solutions as obtained from the stability analysis of such solutions. Whereas the 
soliton solutions are allowed for arbitrary values of the nonlinear parameter, the stability 
condition on the soliton solutions puts restriction on the nonlinearity parameter for which 
stable soliton solutions are allowed [3,4,6,10]. On the other hand, the compacton solutions 
are allowed only within a certain range of the nonlinear parameter (the range is determined 
from the condition of the finite derivative of the compacton solutions at the edges [2,5]) and 
all the allowed compacton solutions (within this allowed range of the nonlinear parameter) 
are stable. Unlike soliton solitions, the stability of the compacton solutions does not put 
any additional constraint on the range of the nonlinear parameter. This result is true even 
for the higher order nonlinear dispersion equations for compactons, whereas for the soliton 
case the higher order linear dispersion term stabilizes the solitons with higher power of 
nonlinearity. It may be noted that we are unable to discuss the question of the stability of 
the compacton solutions as considered by Rosenau et al [1] since their compacton equations 
cannot be derived from a Lagrangian. We however suspect that their compacton solutions 
will also be stable. It would be nice if this can be shown in general. 
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